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The drying dynamics in three dimensional porous media are studied with confocal microscopy. 
We observe abrupt air invasions in size from single particle to hundreds of particles. We show 
that these result from the strong flow from menisci in large pores to menisci in small pores during 
drying. This flow causes air invasions to start in large menisci and subsequently spread throughout 
the entire system. We measure the size and structure of the air invasions and show that they are in 
accord with invasion percolation. By varying the particle size and contact angle we unambiguously 
demonstrate that capillary pressure dominates the drying process. 

PACS numbers: 47.56.+r, 47.57.-s, 47.55. Ca 
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Drying is a ubiquitous natural process; it is of particu- 
lar importance and interest for porous media where there 
are many important applications, including soil drying in 
agriculture, recovery of volatile hydrocarbons from un- 
derground oil reservoirs, spray drying in food and phar- 
maceutical industries, and drying of paint and ceramic 
powders. Drying in porous media is often described the- 
oretically by invasion percolation(IP)[l|,|3,l3|,l4'j and mod- 
eled by immiscible displacement experiments in 2D 5]. 
However, despite its great importance, direct imaging 
of drying in three dimensional (3D) porous media is 
rare, largely due to the difficulty of visualization in 3D 
porous media. Instead, indirect measurements such as 
light scattering and acoustic methods [2|, conductance 
measurement [6], or pressure measurement [T], |8|, |9[ have 
been used. While valuable information is deduced with 
these techniques, many important dynamical features 
of drying, such as the flow pattern, the structure of 
the abrupt bursts and the liquid redistribution, remain 
poorly characterized. To study them, the critical inves- 
tigation by direct imaging is indispensable. 

In this paper, we study drying dynamics in 3D porous 
media by direct imaging with confocal microscopy. We 
observe a strong flow from the menisci in large pores 
to the menisci in small pores during drying. The flow 
causes abrupt air invasions starting from large menisci 
which subsequently spread throughout the entire system. 
The size, structure and dynamics of these invasions are 
explained by 3D invasion percolation coupled with liq- 
uid redistributionj,,^]. The liquid redistribution is con- 
firmed by direct observation. By varying the particle size 
and contact angle we also demonstrate unambiguously 
that capillary pressure dominates the drying process. 

Porous media are prepared by evaporating a high 
concentration colloidal suspension, composed of fluo- 
rescently labeled PMMA particles (density Ppmma = 
1.19g/cm^ and index of refraction ripmma = 1.49) sus- 
pended in decalin(/9dec — .897 g / cm^ , Udec = 1.48). The 
particles and solvent are closely index matched, allowing 



3D visualization through more than one hundred particle 
layers using confocal microscopy. Samples are dried on 
clean glass substrates with the initial volume fractions 
between 30% and 40%. 

The overall drying process is demonstrated by evap- 
orating a drop of suspension(3.8/x^) with particle size 
d = 1.1 ± O.lfim. The drying is divided into two stages: 
The suspension is initially compacted as decalin evapo- 
rates and the colloid volume fraction increases. Due to 
the coffee-ring effect [10[ , the suspension compacts at the 
droplet edge with the randomly-close-packed particles, 
creating a network of pores. Images of horizontal(xy) 
and vertical(xz) cross sections taken at the interface be- 
tween the packed region and free suspension are shown 
in Fig. 1(b). Particles can not move once they are in the 
packed region. Interestingly, the xz section indicates that 
the packed region forms a "ceiling" on the top of the 
droplet due to the lowering of the air- liquid interface from 
evaporation. As solvent evaporates, the dense region at 
the edge grows towards the center and eventually fills the 
entire system (Fig. 1(a), top three images). In the second 
stage, air invades the system. Cracks appear first, indi- 
cating large stresses in the system, then air continues to 
invade the network of pores until all the liquid is evapo- 
rated(bottom three images). It is the second stage which 
is the focus of this study. 

This second stage is characterized by a series of abrupt 
air invasions, or bursts|7l. Isl. lol. Illl. Il2l| . A typical burst is 
identified with confocal microscopy by a bright area sud- 
denly turning black, corresponding to the displacement 
of decalin by air, which is optically mismatched, as il- 
lustrated by the highlighted area in Fig. 1(c). In a single 
step, within 0.07s, an area of about 800 particle area is 
invaded. This abrupt behavior indicates that some large 
stress is driving the process. 

The dominant stress during drying is the capillary 
pressure from the tiny menisci [l3|: 
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FIG. 1: Photographs of drying for d = 1.1 ± O.l/ijn colloidal 
suspension, (a), drying images taken from above at six differ- 
ent times: Os(pict.l), 270s(pict.2), 300s(pict.3), 320s(pict.4), 
365s(pict.5) and 410s(pict.6). We divide the process into 
two stages: compacting stage(top row) and air invading 
stage(bottom row), (b), interface between the packed region 
and free suspension, such as the small box area in pict.l. of 
(a). The left picture is a horizontal slice and the right is a ver- 
tical slice. In the right picture, a "ceiling" of closely packed 
particles near the air-liquid interface is clearly visible, (c), 
image of one large burst. An area of 800 particle size (the 
highlighted area) is invaded within 0.07s. (d), cartoon illus- 
trating the flow from large pore to small pore during drying. 



where a is the surface tension coefficient, 9 is the solvent- 
particle contact angle and r is the radius of curvature 
which is typically r ~ O.IR for randomly close packed 
spheres of radius R. Since crj^f^c — 31mN/m and cos9 ~ 1, 
Eq.(l) gives an enormous pressure: AP ^ 6 x 10^ Pa ^ 
6atm for d = I.I/xto particles! Relative to the ambient 
pressure surrounding the droplet, Pq, the liquid within 
the pores has a very low pressure: Pq — AP. Therefore ev- 
ery meniscus acts like a low pressure pump trying to suck 
liquid from other places. Due to the inhomogeneity of the 
pore sizes, the menisci in small pores can produce lower 
pressure and draws liquid from menisci in large pores, as 
illustrated by the cartoon in Fig. 1(d). We emphasize that 
this flow is quite different from the capillary flow in the 
coffee-ring problem, wherein the faster evaporation rate 
at the droplet edge leads to bulk flow outwards from the 
center to the edge. This flow would move the air-liquid 
interface rapidly through the large pores, creating bursts. 



The bursting process will terminate either when all the 
pores are small enough to balance the capillary pressure 
or when the displaced liquid flows to the nearby menisci 
and reduces the capillary pressure. Previous research in 
2D systems proposed that the later mechanism results in 
an exponential cut-off of burst sizej?], |8|. 




FIG. 2: Observation of strong flow from large pores to small 
pores, (a), d — 2/im sample with d = 0.22/im tracer parti- 
cles. The tracer particles move randomly before the invading 
stage, (b), tracer particle trajectories show a strong flow dur- 
ing invading stage, (c), one sample with small pore and large 
pore domains. The darker domain on the left is composed 
by l.l/im particles and the brighter domain on the right is 
by 2fj,m particles, (d), x-y scans of drying process of the 
two-domain sample. Large-pore domain is invaded first. Sur- 
prisingly, the invasion is restricted in the large-pore domain 
until it has been completely invaded, (e), x-z sections at a 
position similar to the dotted line in (d). 



To visualize this flow, we add a small quantity of tiny 
tracer particles((i — 0.22/im) to a suspension of larger 
spheres with d = 2/im,, as shown in Fig. 2 (a). The tracer 
particles are small enough to move freely through the 
pores. During the compaction stage, the tracer particles 
move randomly J^]. However, upon entering the invasion 
stage, a strong flow appears, as shown in Fig. 2(b). The 
flow velocity and direction vary with time and location, 
with the magnitude measured from 5 to 50^m/s. The 
flow is not from center towards edge, thus we exclude 
the possibility of the "coffee-ring" effect. This provides 
direct evidence of the flow in Fig. 1(d). 

To conflrm that the flow is indeed from large to small 
menisci, we prepare samples comprised of two distinct do- 
mains of differing particle diameters (d=l.l and 2 um), il- 
lustrated in Fig. 2(c). Various times in the drying process 
at the boundary of the two domains are highlighted by 
the xy sections shown in Fig. 2(d). The large-pore domain 
(right) is invaded first, proving that air invasions start in 
large pores. Strikingly, fluid displacements are restricted 
to the large-pore domain until it has been completely 



invaded. Since evaporation occurs everywhere, the only 
way to halt air invasions at the boundary between the do- 
mains is to cause a macroscopic flow from the large-pore 
to the small-pore domain to compensate for evaporation. 
The large-pore domain effectively serves as a liquid reser- 
voir which prevents drying of the small-pore domain until 
the reservoir is completely exhausted. A similar process 
in the xz direction is shown in Fig.2(e). This phenomenon 
is further proof confirming the flow pattern from large to 
small menisci shown in Fig. 1(d). Moreover, it provides 
a means to precisely control the order of drying. It also 
illustrates a method to spontaneously concentrate mate- 
rial to small-pore regions through drying. 
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FIG. 3: Size, structure and dynamics of bursts. Main panel 
of (a), the probability, P, of finding a burst of size s or 
A. Both the 2D (•) and 3D (o) size distributions are plot- 
ted. The 3D size, s, is computed from the 2D size. A, 
with the relationship s ^ A f' . The two fitting func- 
tions are: P{A) ~ A~°'Exp{-A/A*), with a = 1.6 ± 0.1, 
A* = 850 ± 150/im^ and P{s) ~ s'^' Expi-s/s"), with 
r' = 1.5 ± 0.1, s* = 6500 ± 1000. The smallest area in plot is 
limited by the image resolution. Lower inset, zoom-in picture 
of freshly invaded region. Only 15% pore space is filled by 
air. Upper inset, global pattern by many invasions. The field 
of view is 345/im x 345/im. (b), liquid redistribution due to 
a burst. The 1st image is before the burst and The 2nd is af- 
ter. An invasion occurs at the lower arrow spot, resulting in 
an obvious meniscus re-adjustment at the upper arrow spot. 
The 3rd picture traces the air-liquid interface of the first two 
images. The dotted curve is before the burst and the solid one 
is after. It shows menisci re-adjustments at several positions. 



The flow evacuates the liquid from large pores and 



leads to abrupt bursts. The sequence of bursts has 
been studied indirectly through the associated pressure 
changes in similar systems [7|,|8|; however, there has been 
no structural characterization of these events. We di- 
rectly visualize these bursts by horizontally scanning a 
fixed area over time. The horizontal scanning requires 
drying to proceed horizontally as well. We achieve this by 
drying the sample between two parallel substrates. We 
coat both substrates with a layer of PMMA to match the 
wetting property with the bulk particles. Since the in- 
vaded air strongly scatters the light, we restrict our mea- 
surements to regions near the bottom substrate, where 
the measurement is most reliable. Images are acquired 
at 2.5 frames per second, and the burst areas are mea- 
sured by subtracting two adjacent frames (see fig 1(c)). 
More than 600 bursts in a single drying process are mea- 
sured to obtain good statistics. We find that the area 
of the bursts ranges from a single particle to hundreds 
of particles. The probability of a burst, P, as a func- 
tion of its area. A, exhibits a power-law distribution for 
small areas but rapidly decreases for large areas as shown 
by the log-log plot(») in the main panel of Fig. 3(a). The 
power-law distribution is in accord with the invasion per- 
colation and the deviation at large areas reflects the fact 
that, during the bursts, liquid does not have time to evap- 
orate but instead gets redistributed among the menisci. 
This redistribution decreases the capillary pressure and 
prevents the bursts from growing too large. The func- 
tional form P{A) '- ^~"ii;xp(-A/^*)0,U, accounts for 
both the power-law behavior of invasion percolation and 
the exponential decay due to redistribution of the liquid 
and fits the data well as shown by the solid line in Fig. 3. 
We obtain a = 1.6 ± 0.1 and A* = 850 ± 150^m2. 

The exponent differs from previous work which found 
a = 1.32 [71, |8[. However, this was for a 2D network of 
pores whereas our system is roughly one hundred par- 
ticle layers in thickness and thus more nearly 3D. We 
examine the freshly invaded regions [l5J at single parti- 
cle level and find a complex structure, as shown in the 
lower inset of Fig. 3(a). Instead of being completely filled 
by air, only 15 ~ 25% of the invaded region is filled by 
air. This value agrees with 3D bond percolation thresh- 
olds, 12 ~ 25% [165 . Scans at different heights also show 
that the invaded region is a 3D connected network of 
air. These results suggest that each burst is a 3D inva- 
sion percolation event. Therefore we compute the true 
3D size of a burst, s, from the 2D area. A: s ~ A-^f^^, 
with Df ~ 2.5 the fractal dimension of 3D percolation. 
We plot P{s) in the main panel as open symbols(o) and 
fit the data with: P{s) ^ s~'^ Exp{—s/s*), and obtain 
t' = 1.5±0.1, s* = 6500±1000. Theoretical work showed 
that the power r' is related to percolation exponents [1 71]: 
r'Df — Df + De — 1^^^- This relation predicts r' = 1.55 
in 3D and 1.31 in 2D. Our data are in excellent accord 
with the 3D result. 

The sequence of bursts quickly spreads throughout the 



field of view and forms a macroscopic pattern, as shown 
in the picture taken at the moment when invasions first 
percolate the field of view(345/i7TT. x 345/iTO) in the upper 
inset of Fig. 3 (a). The pattern is fractal- like and similar 
to previous experiment [3| • However, from the magnified 
image in lower inset, we determine that the black region 
is not completely filled by air, but instead formes 3D 
percolated networks of air. 

We can also directly determine the origin for the expo- 
nential cut-off: the liquid redistribution. In two frames 
separated by 0.07s, we catch a single burst and the subse- 
quent liquid redistribution. The two frames are shown in 
the two pictures of Fig. 3(b) on the left. The first picture 
is before the burst while the second is after. The burst oc- 
curs at the spot indicated by the lower arrow and an obvi- 
ous meniscus re-adjustment occurs at the spot indicated 
by the upper arrow. More careful comparison is done 
in the third picture of Fig. 3(b), where we trace the air- 
liquid interfaces of the first two pictures. The dotted line 
is the interface before the burst while the solid line is af- 
ter. There are meniscus re-adjustments in several places. 
These re-adjustments reduce the capillary pressure and 
lead to the exponential cut-off, s* . This observation pro- 
vides unambiguous evidence for the redistribution argu- 
ment proposed by previous indirect measurements [7|,|8|. 

We also vary both the pore size and the contact angle, 
to further validate the prediction of Eq.(l). The effect of 
pore size is demonstrated by comparing two completely 
dried samples with particles of d = 0.22^m{leit) and 
d — 2/ir7i(right) in Fig. 4(a). The beautiful crack pattern 
in the d — 0.22/im sample indicates large stresses dur- 
ing drying, while the absence of cracks in the d = 2fim 
sample suggests that the stresses are much lower; This 
is consistent with Eq.(l). We vary the contact angle of 
the invading fluid by studying the re-hydration of a dried 
sample. We see a marked difference in re-hydration: the 
liquid-air interface sweeps the whole field of view in a 
smooth and continuous manner, in sharp contrast with 
the fractal pattern of drying in Fig. 3(a). This demon- 
strates the important effect of 9, as predicted in Eq.(l). 
This result is in great agreement with a previous simula- 
tion, which explained the smooth interface by the coop- 
erative activity of neighboring menisci | ISl. 1 19l|. 

By using direct imaging, we show that drying of a 3D 
porous media can be well described by invasion perco- 
lation provided the effects of liquid redistribution are 
included; which leads to an exponential cut-off of the 
power-law distribution of abrupt invasions. The liquid 
redistribution is directly observed, verifying the impor- 
tant feature that had heretofore only been deduced by 
indirect methods. Our study reveals the flow pattern of 
drying which may ultimately afford a means to control 
drying in porous media. A deep understanding of drying 
may provide a means to help eliminate cracking or other 
undesirable effects. 

This work was supported by ICI with partial sup- 
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FIG. 4: Variation of particle size and contact angle, (a), 
completely dried samples of small (d — G.22^m, left) and 
large (d = 2^,m, right) particles. There are many cracks for 
0.22^m sample and no cracks for 2^m one. (b), re-hydration 
of a dried sample. We see a smooth and continuous interface, 
in sharp contrast to the fractal pattern of drying. 
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